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Abstract 

We study the statistical properties of overdamped particles driven by two cross-correlated multi- 
plicative Gaussian white noises in a time-dependent environment. Using the Langevin and Fokker- 
Planck approaches, we derive the exact probability distribution function for the particle positions, 
calculate its moments and find their corresponding long-time, asymptotic behaviors. The generally 
anomalous difi'usive regimes of the particles are classified, and their dependence on the friction 
coefficient and the characteristics of the noises is analyzed in detail. The asymptotic predictions 
are confirmed by exact solutions for two examples. 

PACS numbers: 05.40.-a, 05.10.Gg, 02.50.-r 



1 



I. INTRODUCTION 



Diffusive behavior is intrinsic to many physical, chemical, biological, economical and 
other systems ll^Sy]- One of the most rigorous approaches for its description that takes 
into account the dynamical origin of diffusive motion is based on the Langevin equation, 
i.e., the stochastic equation of motion. For a large variety of systems that are stochastically 
equivalent to free Brownian particles, this equation has the simplest form, namely a constant 
friction coefficient and additive Gaussian white noise. Such systems exhibit normal diffusion, 
i.e., the long-time asymptotics of the dispersion of the state parameter are proportional to 
time. 



However, plenty of systems exhibit anomalous diffusion (for a review, see Refs. 



In general, there exist several approaches to describe this phenomenon. A first one is based 



on fractional calculus and involves fractional diffusion equations 



Planck equations 
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14l | or fractional Langevin equations 



fractional Fokker- 
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17l |. Yet another 



other relies on generalized ordinary Langevin equations. The latter is very attractive for 
the treatment of diffusive behavior, and it is especially informative if those equations can 



be solved exactly (18|, |19| . This approach has been employed to study anomalous diffusion 
for a variety of systems, including systems described by the generalized Langevin equation 
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23L |2y| , free undamped 



with a friction memory kernel j2 

damped particles js^ driven by additive noise, overdamped particles driven by one 
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and 



31( 1 and 



two multiplicative colored noises, and overdamped particles with time-dependent drift 
driven by additive Gaussian white noise js^. 

The Langevin approach is also an efficient tool to study various effects in systems driven 
by correlated noises. The correlation can exist, for example, if the two noises possess a 



common origm 



35|, as in the case of bistable systems where fluctuations of some model 



parameters are equivalent to the action of correlated additive and multiplicative noises [3, 



Such noises play an important role in the description of stochastic resonance 
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non- 



Q- We expect that 



equilibrium phase transitions and noise-induced transport 

correlated noise will impact significantly the anomalous diffusive behavior. To investigate 
this problem in the general case of a time-dependent environment, we consider the long-time 



2 



behavior of an overdamped particle governed by the Langevin equation 

2 

mm = Y.gMt))m [x{o) = xoi (i.i) 

i=l 

where x{t) is the particle position, A(t) denotes an arbitrary, positive- valued time- dependent 
friction coefficient, and gi{x) are deterministic functions that characterize the state- 
dependent action of Gaussian noises C,i{t). We suppose that all parameters and variables in 
Eq. are dimensionless and that the noises are white with (^j(t)) = and 

mMt'))=2^^At-t')■ (1-2) 

Here (. . .) denotes averaging with respect to the noises ^i(t), An = Ai(> 0) and A22 = 



A2(> 0) are the intensities of the noises ^i(t) and ^2(^)5 respectively, A12 = A21 = rvAiA2, 
r is the parameter characterizing the cross-correlation of the noises, \r\ < 1, and 6{t) is 
the Dirac delta-function. For generality, we assume that the noises ^i(t) are external and 
so the time-dependent friction and stochastic forces in Eq. p.lj) are not related to the 
fluctuation-dissipation theorem. Hence, if internal noises are negligible in comparison with 
external ones, this equation can be used, e.g., for the description of Brownian particles in a 
liquid with a time-dependent temperature. It is important to note, however, that Eq. 
has a physical meaning for internal noises as well. Indeed, since according to the Stokes 
formula the friction coefficient depends on the particle radius, this equation describes in 
the overdamped regime the Brownian dynamics of spherical particles with time-dependent 
radii. Therefore, Eq. (jl.ip with internal noises can be used, e.g., for the study of stochastic 
dynamics of nucleating and reacting particles, see in this context also related works and the 



given applications therein which considers instead a fluctuating friction |4J, 44 1. Equation 
(II. 1|) does not include inertial effects since in most cases they are important only for a very 
short initial time interval and can be safely neglected for larger time scales. 

The paper is organized as follows. In Sec. II, using the Langevin approach, we derive the 
time-dependent probability distribution function (PDF) of x{t) and calculate its moments. 
In Sec. HI, we flnd the long-time asymptotic behavior for moments and classify the diffusive 
regimes of the particles. We obtain some exact solutions in Sec. IV. In the same section we 
compare the exact and asymptotic results. We summarize our results in Sec. IV. Finally, we 
derive the PDF by the Fokker-Planck method in the Appendix. 
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II. PROBABILITY DISTRIBUTION FUNCTION AND ITS MOMENTS 



According to the results of Refs. [3J, |35|, |45[ , the two- noise Langevin equation ()1.H) is 



statistically e quiv alent to a one-noise Langevin equation. Its form depends on how Eq. 
is interpreted j4o| . In the case of the Stratonovich interpretation j46,|47| we obtain 

X{t)x{t) = G{x{t))C{t) [x{0) = xo], (2.1) 

where 



2 2 



J2J2^^M^)9A^) (2-2) 
\ i=i j=i 



G{x) 

and ({t) is Gaussian white noise with zero mean and correlation function 

{cmt')) = mt-t'). (2.3) 

In general, the function G{x) will be nonnegative. Here, we consider the slightly more 
stringent case that G{x) > for all x, and consequently x(t) G {—oo, oo). 

In order to find the PDF of x{t), we use an explicit expression for x{t) that follows 
from Eq. (j2.1|) . (An alternative derivation of the PDF via the solution of the corresponding 



47 



Fokker-Planck equation is presented in the Appendix.) The Wang and Zakai theorem 
states that the conventional rules of calculus are applicable to Langevin equations if they are 
interpreted as Stratonovich stochastic equations. Accordingly we solve for x(t) by separating 
the variables and integrating Eq. p.lj) to obtain 

K*) dx' fH{t')dt' 



XO 



G{x') Jo A(t') 



(2.4) 



If we define 

^{x)-^>{x,)= / -— (2.5) 



G{x') 



and 



then Eq. ()2.4|) can be rewritten in the form 

^{x{t)) = w{t) + Uo (2.7) 

[uq = \E'(xo)]. According to the definition ()2.5p . the function u = \E'(x) is continuous and 
monotonically increasing; consequently, we obtain minw = \1'(— oo), maxu = \E'(oo). On 



the other hand, the range of the function w{t) is over all real- valued numbers. Therefore, 
Eq. ()2.7|) is applicable at all times only if \E'(— oo) = — cxd and \E'(oo) = oo. We assume 
that these conditions are met and take into account that the inverse function x = \E'^^(m) is 
single-valued, i.e. it is also continuous and monotonically increasing, to obtain the unique 
solution of Eq. (j2.H) in the form 

x{t) = ^-\w{t) + uo). (2.8) 

The relation ()2.8j) shows that the statistical properties of x{t) are defined by the Gaussian 
process w{t), which is fully characterized by its zero mean and the correlation function 

{w{t)w{t')) = 2 (2.9) 

Let Px{x,t) and Pw{w,t) be the univariate PDFs that x{t) = x and w{t) = w, respectively. 
Since a one-to-one correspondence exists between x{t) and w(t), the condition Px{x, t)\dx\ = 
Pwiw,t)\dw\ must hold. Taking into account that Pw{w^t) = [-\/27r a"(t)]^^ exp[— tf^/2cr^(t)], 
where 

a'{t) = 2f'-^ (2.10) 



/o A2(rj 

is the dispersion of w{t), and using the relation that dx/dw = G{x), which follows from 
Eqs. fj2.5|) and ()2.7|) . we obtain for the desired PDF: 



PxixA) = —= exp 

^ ' V2^G{x)a{t) 



[^{x)-uoY 



(2.11) 



2^2 (t) 

This distribution function is properly normalized, namely, Px{x,t)dx = 1, and its 
moments, defined as 

/oo 
dxx'^Pxix.t) (2.12) 
-oo 

(n = 1, 2, . . .), contain valuable information about the stochastic system dynamics. Substi- 
tuting Eq. ()2.11|) into Eq. ()2.12|) and introducing the new variable y = {"^{x) — uq) j a{t), we 
reduce Eq. fl2.12j) into the convenient form 

(x"(t)) = ^/ dyh>-\ya{t)+u,)\ e"^ (2.13) 



/27r 

Depending on the time-dependence of the friction coefficient A(t), the particles can be 
localized or delocalized as t — > cxo. If a{oo) < oo, i.e., if A(t) grows fast enough, then the 
stationary PDF Px{x,oo) exists, and the particles are localized. Brownian particles in a 
freezing liquid serve as an example. On the contrary, ii a(t) oo as t — ^> oo, then the 
particles are delocalized. We focus on the latter situation in the following. 



III. ASYMPTOTIC BEHAVIOR FOR THE MOMENTS 



To study the long-time behavior of x{t), including various diffusive regimes, we derive 
the asymptotic behavior of (a;"(t)) for t ^ oo. This is a rather complex problem because 
the integrand in (|2.13j) contains the inverse function \I'~^('u) of the integral function 
Since the dispersion cr^(t) tends to infinity in the long-time limit, the inverse function at 
y > and y < can be replaced by its asymptotics at n ^ oo and u — oo, respectively. 
To find these asymptotic behaviors, we assume that 



G{x) jAplxl""" (x-^poo). 



(3.1) 



where p = + or — , Ap is the effective white noise intensity at x = p oo, and ap < 1 to ensure 
that \I'(±oo) = ±oo. Then for x ^ poo, Eq. ()2.5|1 yields 



In—, Op = 1, 



I ^ I 1 Oip 

1- an 



(3.2) 



an < 1. 



[The normalizing parameter Op can be found if the explicit form of the function \l/(x) is 
known (see Sec. IV).] Finally, from Eq. ()3.2p we obtain 



\E' ""^(m) ~ papexp(w Ap |m|) (w-^poo) 



(3.3) 



for Op = 1, and 



"^'^u) p[Ja'p{1 - ap)\u\Y-°'p (w-^poo) 



(3.4) 



for ctp < 1. Next we consider separately these two cases. 



1. Up = 1 



In this case, we replace the inverse function ^ ^{ya{t) + uq) by its asymptotic behavior, 
i.e., a_|_ exp[-\/A+(y(T(t) + uq)] at y > and — a_ exp[— A/A_(y(T(t) + uq)] at y < 0, and 
Eq. (jTT^ thus yields 



1 r°° I — 

(x"(t)) ~ -= E(P«p)" / dy exp(n Ja; {ya{t) 



+ PMo) -1/V2) 



(3.5) 
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as t — > oo. Using the standard integral formula |4t 

fdye-^y-y'/' = y|e-^/^erfc(^), (3.6) 

where erfc (2;) = (2/y^) J^dte~^^ is the complementary error function, and taking into 
account the relation erfc(— 00) = 2, we reduce the asymptotic formula ()3.5p to the form 

(x"(t)) ~ E(P«p)" exp{pn^,Uo + /\pa\t) /2). (3.7) 
p 

It is important to note that, since ()3.3|) represents the leading term of the asymptotic 
expansion of \I^^^(u), we should keep only the largest term in the right-hand side of ()3.7p . 
In particular, if A-(. > A_ then 

(x"(t)) ~ al e^vin^X^Uo + /\+a'^{t) /2) (3.8) 

{t 00), i.e., the particles tend to plus infinity, and if A+ < A_ then 

(x"(t)) ~ (-l)"a!! e^^{-n^XluQ + n"" A^a^t) /2) (3.9) 

{t 00), i.e., the particles tend to minus infinity. We emphasize that such a behavior of 
x{t) results as a direct consequence of the multiplicative nature of the noise. 

If A+ = A_ = A then a_|_ = a_ = a. Both terms in the right-hand side of (|3.7|) have the 
same order, and the asymptotic relation (j3.7j) is reduced to 

(x"(t)) ~ 2a" cosh(riVAuo) exp(n^Aa2(t)/2) (3.10) 

(t — > 00) for even n and to 

(a;"(t)) ~ 2a'^sinh(nVAMo) ex.^{n^ /^a"^ {t) / 2) (3.11) 

{t —>■ 00) for odd n. The last relation shows that, in contrast to the previous discussed cases, 
no systematic growth of {x{t)) occurs for uq = 0, i.e., xq = \I^~^(0). In other words, in this 
case the particle position x(t) exhibits purely diffusive behavior that can be characterized 
by the dispersion cr^(t) = {x^(t)) — {x(t))'^. Since al{t) ~ 2a^e^^'^^(*^ as t — > 00, different 
long-time regimes of diffusion can exist. 

We describe these regimes for the case that X{t) ~ Zt^ {t —>■ 00), where / is a positive 
parameter and /? < 1/2 (the last condition guarantees that a (00) = 00). Using Eq. (|2.1U|) . 
for /? < 1/2 we obtain 



{t oo). Accordingly, we call the diffusive behavior a stretched exponential one if < 
P < 1/2, exponential if /3 = 0, and a compressed exponential if /3 < 0. If /3 = 1/2, then 
cr^(t) ~ (2//^) ln(t/t) [t is some characteristic time scale] and so 

cr^(t)~2a2(t/t)t (3.13) 

as t ^ oo. This asymptotic formula demonstrates a truly remarkable feature of this physical 
system, namely: the character of its diffusive regime depends on the effective white noise 
intensity A which in the general case depends on the correlation coefficient r. Following the 
conventional terminology, we say that the system exhibits subdiffusion if A < /^/4, normal 
diffusion if A = P/A, and superdiffusion if A > P/A. We expect that these diffusion regimes 
can exist in the systems of nucleating particles whose radii grow as t^^"^. 



2. Op <1 

In this case, using the standard integral formula j^] 

f <^y^/'•e-»v^ = 2'^'^(^i±i) (3,i4) 

> —1, r(2;) = dx x^~^e~^ is the gamma function] and the asymptotic relation (j3.4j) . 
Eq. fll^ yields 

(:^"W)~^E(pirr(W2) 

2V7r 



x[^pil- a,)a{t)]^^-' (3.15) 



(t oo), where rjp = 1 + n/{l — Op). We emphasize that in this asymptotic relation only 
the dominant term on the right side should be kept. Specifically, if a+ > a_ then 

(^"W) - ^^^^[s/2A;{1 -a^)a{t)r~\ (3.16) 

and if a+ < a_ then 



{x-{t)) ~ {-ir^^^^^[^{l~a.)a{t)r-\ (3.17) 
Finally, if = a {r]p = rf) then both terms have the same order and so 

V{r}l2) nz± 

(x"(t)) ~ ^f^[A+^ + (-1)"A J ] 

x[v^(l -a)a(t)]"-^ (3.18) 
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These relations show that for A_|_ 7^ A_ all moments of x{t) diverge as t — 00. 

Next, we assume that ap = a and Ap = A. In this case, using the asymptotic relation 
(Q, Eq. ^rn^ for t cx) is reduced to 

~ IM^[V2A (1 - a)a{t)r' (3.19) 



(a;-(t)) ^ r f ^ ) [v^ (1 - a)]"- V^-2(t) (3.20) 



if n is even, and to 



\/2^ V 2 

if n is odd. The last relation and the exact formula (j2.13p confirm the fact that {x{t)) = 
for uq = 0, i.e., under these conditions the particles display diffusive behavior. Since in this 

case 

, r(l/(l - «) - 1/2), ^ .^2 , , 

alit) ~ ^ ^[v^(l - a)cr{t)]— (3.21) 



as t 00, its character depends on the exponent a and on the asymptotic behavior of 
a{t). Specifically, if \{t) oc {t ^ 00, p < 1/2) and so a^{t) oc t^-^/? p < and 
(7^(i(:) oc ln(t/t) for P = 1/2, then the particles exhibit super diffusion if a > 2/3 and P < 1/2, 
normal diffusion if a = 2/5 and /? < 1/2, subdiffusion if a < 2/5 and /? < 1/2, and logarithmic 
diffusion [al{t) oc lnT^(t/t)] if /5 = 1/2. 

Thus, multiplicativity of noises and time dependency of the friction coefficient can give 
rise to anomalous behavior of the system. Specifically, if A = const then the conditions 
G{x)\x^±oo 00 and G{x)\x^±oo are responsible for the fast and slow diffusion, 
respectively. On the contrary, if G{x) = const (white noises are additive) then diffusion 
is fast if X{t)\t^oo and is slow if A(t)|(^oo 00. Remarkable, the presence of both 
these factors can lead to new stochastic phenomena, like the dependence of the character 
of anomalous diffusion on the effective white noise intensity. We note also that the nature 
of anomalous diffusion in this system, due the features of the stochastic and friction forces 
mentioned above, is quite different from that observed for random walks with long-tail 
jump-length and/or waiting-time distributions. 



IV. SPECIFIC EXAMPLES 



To confirm our asymptotic results, we calculate explicit the expressions for the moments 
of x{t) for two particular cases. 
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A. Two additive noises 



As a first example, we consider tlie case that both white noises are additive, i.e., gi{x) = 
g2{x) = 1. Then, G\x) = Ai + A2 + 2ry^A^ = G^, ^(x) = x/G, ^-^(u) = Gu, 
uq = xq/G, and Eq. ()2.13|1 is reduced to 

(x"(t)) = ^ / dy [yGait) + xoT e'^'/^ (4.1) 

V ZTT J -00 

Using in 1)4.11) the binomial formula and taking into account the relation ()3.14|) . we obtain 
after some calculations 



[n/2] . 1 



V2G 



a{t) 



2k 



(4.2) 



where [n/2] is the integer part of n/2 and = n\/{n — m)!m! is the binomial coefficient. 

According to this formula, the leading term of the long-time asymptotic expansion of 
{x^{t)) is given by the term in the sum that corresponds to = [n/2], i.e., 

(x"(t)) r^^T(^][V2Ga{t)r (4.3) 



> V 2 
for even n and 

(x"(t))~^nr(^)[v^Ga(t)]"-i (4.4) 

for odd n. Since A = G^, a = and rj = n + 1, the same asymptotic behavior follows also 
from the formulas ()3.19p and ()3.2U|) . respectively. 

B. Additive and multiplicative noises 

As a second example, we consider a linear Langevin equation where one noise is multi- 
plicative, gi{x) = X, and the other one is additive, i.e., g2{x) = 1. According to the definition 
()2.2|1 . in this case 



G(x) = A/Ai(x2 + 2rz/x + z/2) (4.5) 



(z/ = y A2/A1 ) and for |r| < 1 we obtain 

1 f X + ru \ 

W[x) = ^= arcsinh — (4.6) 
V Ai Vz/y 1 — r^/ 

and 

$-1(m) = uVl - r2sinh(y^M) - ru. (4.7) 
10 



Then, using Eqs. fj2.13|) and ()4.7p . the relation 



sinh(A/Ai Mo) 



V 



(4.8) 



that follows from Eq. ()4.7|) . and the standard integral formula 



dy e 



-hy^+cy 



cosh(p?/). 
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X 



sinh(cp/26) 
cosh(cp/26). 



(4.9) 



(6 > 0), we find the first moment 



(x(t)) = (a;o + rv) exp(Aio-^(t)/2) - rz/, 



(4.10) 



which shows that x{t) exhibits purely diffusive behavior only if r = —xq/v, and the second 
one 

{x\t)) = [(xo + ruf + u\l- r^)/2\ exp[2Ai(T'(t)] 
- 2rz/(xo + rz/) exp(Ai(T^(t) /2) 

+ z/2(3r2 - l)/2 . (4.11) 



For t — >^ cxD, these exact results readily yield 



{x{t)) ~ (xo + rz/) exp(Ai(T2(t)/2) 



(4.12) 



and 



(x=^(t)) ~ [(xo + rz/)2 + z/2(l - r2)/2] exp[2Aia2(t)]. (4.13) 

On the other hand, the same long-time representations for {x(t)) and {x^{t)) follow as well 
from the more general asymptotic formulas ()3.1()|1 and 1)3.111) . To verify this, we compare 
()3.1|) with the asymptotic relation G{x) ~ v^A7 |x|, which follows from Eq. ()4.5|) as |x| ^ cxd; 
this yields Ap = A = Ai and ap = a = 1. Next, since arcsinh^; ~ pln2|z| as z ^ poo, 
from Eq. ()4.6|) we find the asymptotic formula 

21x1 



*(x) 



_P_ 
/A7 



In 



(4.14) 
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(x poo), which by comparison with ()3.2p yields ttp = a = v\J\ — r'^ /2. Finally, substitut- 
ing the values for A and a into the formulas (jH.llj) [with n = 1] and [with n = 2] and 
using the relation ()4.8j) . we indeed recover the asymptotic behavior ()4.12j) and ()4.1Hj) . 

Thus, our exact results corroborate the general asymptotic formulas obtained in the 
previous section. 

V. CONCLUSIONS 

We have studied analytically the statistical properties of a special class of exactly solv- 
able stochastic models represented by an overdamped particle that is driven by two cross- 
correlated multiplicative Gaussian white noises in a time-dependent environment. As a first 
step, we have reduced the initial two-noise Langevin equation that describes the particle 
dynamics in the Stratonovich sense to a stochastically equivalent one-noise Langevin equa- 
tion. Then, solving the latter equation and the corresponding Fokker-Planck equation, we 
have derived the exact probability distribution function for the particle positions. 

To study the long-time behavior of the particle dynamics, we have calculated the asymp- 
totic behavior for the moments of the distribution function. Their analysis results in a rich 
behavior of differing regimes of anomalous particle diffusion; namely we find regimes with 
normal diffusion, subdiffusion, super diffusion, exponential diffusion, stretched exponential 
diffusion, compressed exponential diffusion, and, as well, logarithmic diffusion. We have 
established conditions for these diffusion behaviors to occur and in some special cases we 
were able to derive exact formulas for the moments of order 2 and higher that confirmed the 
general asymptotic results. Also, we have discovered the truly remarkable feature for sys- 
tems described by a special class of linear Langevin equations that their diffusive behavior 
becomes determined by the effective white noise intensity. 
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APPENDIX: ALTERNATIVE DERIVATION OF THE PDF (fCT^ 



To derive the PDF ()2.1H) by the Fokker- Planck method, we assume that each Gaussian 
white noise ^j(t) is characterized by its own parameter 7^ (0 < 7j < 1) that determines the 
points of time at which gi{x{t)) is evaluated in the corresponding integral sum. Then the 
Fokker-Planck equation associated with the Langevin equation takes the form 

\\t)—P^{x, t) = -—h{x)P^{x, t) + —d{x)P^{x, t) (A.l) 



[Px{x, 0) = 5{x — Xq)], where 



h{x) = 2j2jZl^^id{^)gA^) (A.2) 
i=l i=l 



is the drift coefficient, 

2 2 

d{x) = Y.Y. ^ijgi{x)gj{x), (A.3) 
i=i j=i 

is the diffusion coefficient, and the prime denotes the derivative with respect to the argument 
of the function. 

In order to solve Eq. ()A.1|) . we introduce a new variable u = U{x), where the function U{x) 
remains to be defined. If the functions U{x) and U^^{u) are single- valued and U'{x) > 0, 
then 

Pxix,t) = P^{u,t)U'ix). (A.4) 

According to Eq. ()A.1|) . the PDF Pu{u,t) of the random process u(t) = U{x(t)) satisfies 
the Fokker-Planck equation whose drift and diffusion coefficients (as functions of the old 
variable x) are given by 

h{x) = h{x)U'{x) + d{x)U"{x), d{x) = d{x)[U'{x)]^. (A.5) 



If we chose d{x) = 1, then U'{x) = l/yd(x) and 

~ , , 2h(x) — d'ix) , . , 

h{x) = ^ ' ^ -. (A.6) 

2,Jd{x) 

Remarkably, if 7^ = 1/2, i.e., if the Stratonovich interpretation of the Langevin equation 
fll.l|) is used, then h{x) = and the Fokker-Planck equation for Pu{u,t) takes the simplest 
form 

X\t)-Puiu,t) = —P.,iu,t). (A.7) 
13 



Solving it, for example, by the generating function method and using the initial condition 



P(m, 0) = 5{x — xo)\/d{x) [recall that x = U ^{u)], we find 



Pu{u,t) 



■ exp 



[u - Up) 
2a2(t) 



(A.8) 



/2na{t) 

Finally, substituting this result into Eq. ()A.4jl and taking into account that U{x) = 
d{x) = G'^{x) and -u = \E'(x), we again obtain Eq. ()2.1H) . 



[1] N. G. Van Kampen, Stochastic Processes in Physics and Chemistry (North-Holland, Amster- 
dam, 1992). 

[2] W. Horsthemke and R. Lefever, Noise-Induced Transitions (Springer- Verlag, Berlin, 1984). 
[3] C. W. Gardiner, Handbook of Stochastic Methods, 2nd ed. (Springer- Verlag, Berlin, 1990). 
[4] J. P. Bouchaud and A. Georges, Phys. Rep. 195, 127 (1990). 

[5] D. ben-Avraham and S. Havlin, Diffusion and Reactions in Fractals and Disordered Systems 

(Cambridge University Press, Cambridge, 2000). 
[6] G. Zaslavsky Phys. Rep. 371, 461 (2002). 

[7] R. Metzler and J. Klafter, J. Phys. A: Math. Gen. 37 R161 (2004). 
[8] W. R. Schneider and W. Wyss, J. Math. Phys. 30, 134 (1989). 

[9] R. Hilfer, Fractals 3, 211 (1995); R. Hilfer, Fractional calculus in physics, (World Scientific, 
Singapore, 2000). 
[10] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000). 
[11] H. C. Fogedby Phys. Rev. E 50, 1657 (1994). 
[12] A. Compte, Phys. Rev. E 53, 4191 (1996). 

[13] R. Metzler, E. Barkai, and J. Klafter, Phys. Rev. Lett. 82, 3563 (1999). 
[14] E. Barkai, Phys. Rev. E 63, 046118 (2001). 

[15] V. Kobelev and E. Romanov, Prog. Theor. Phys. Suppl. 139, 470 (2000). 
[16] E. Lutz, Phys. Rev. E 64, 051106 (2001). 

[17] A. Bazzani, G. Bassi, and G. Tmxhetti, Physica A 324, 530 (2003). 
[18] P. Hanggi, Phys. Lett. A 83, 196 (1981). 

[19] B. O'Shaughnessy and I. Procaccia, Phys. Rev. Lett. 54, 455 (1985). 

[20] R. Murahdhar, D. Ramkrishna, H. Nakanishi, and D. Jacobs, Physica A 167, 539 (1990). 

14 



[21] J. M. Porra, K. G. Wang, and J. Masoliver, Phys. Rev. E 53, 5872 (1996). 

[22] K. G. Wang and M. Tokuyama, Physica A 265, 341 (1999). 

[23] R. Morgado, F. A. Oliveira, G. G. Batrouni, and A. Hansen, Phys. Rev. Lett. 89 100601 

(2002). 

[24] J. D. Bao and Y. Z. Zhuo, Phys. Rev. E 71, 010102(R) (2005). 

[25] J. D. Bao, P. Hanggi, and Y. Z. Zhuo, Phys. Rev. E 72, 061107 (2005). 

[26] A. V. Mokshin, R. M. Yulmetyev, and P. Hanggi, Phys. Rev. Lett. 95, 200601 (2005). 

[27] J. Masohver, Phys. Rev. A 45, 706 (1992). 

[28] J. Hcinrichs, Phys. Rev. E 47, 3007 (1993). 

[29] J. Masohver and K. G. Wang, Phys. Rev. E 51, 2987 (1995). 

[30] S. I. Denisov and W. Horsthemke, Phys. Rev. E 62, 7729 (2000). 

[31] S. I. Denisov and W. Horsthemke, Phys. Rev. E 65, 031105 (2002). 

[32] A. N. Vitrenko, Physica A 359, 65 (2006). 

[33] F. Lillo and R. N. Mantegna, Phys. Rev. E 61, R4675 (2000). 

[34] A. Fuhnski and T. Telejko, Phys. Lett. A 152, 11 (1991). 

[35] A. J. R. Madureira, P. Hanggi, and H. S. Wio, Phys. Lett. A 217, 248 (1996). 

[36] L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni, Rev. Mod. Phys. 70, 223 (1998). 

[37] P. Hanggi, ChcmPhysChcm 3, 285 (2002). 

[38] J. H. Li and P. Hanggi, Phys. Rev. E 64, 011106 (2001). 

[39] Y. Jia and J. R. Li, Phys. Rev. Lett. 78, 994 (1997). 

[40] S. I. Denisov, A. N. Vitrenko, and W. Horsthemke, Phys. Rev. E 68, 046132 (2003). 

[41] L. Cao and D. J. Wu, Phys. Rev. E 62, 7478 (2000). 

[42] J. H. Li, J. Luczka, and P. Hanggi, Phys. Rev. E 64, 011113 (2001). 

[43] J. Luczka, P. Talkner, and P. Hanggi, Physica A 278, 18 (1999). 

[44] J. Luczka, P. Hanggi, and A. Gadomski, Phys. Rev. E 51, 5762 (1995). 

[45] D. J. Wu, L. Cao, and S. Z. Ke, Phys. Rev. E 50, 2496 (1994). 

[46] R. L. Stratonovich, SIAM J. Control 4, 362 (1966). 

[47] P. Hanggi and H. Thomas, Phys. Rep. 88, 207 (1982); see Sec. 2.4. 

[48] E. Wong and M. Zakai, Ann. Math. Stat. 36, 1560 (1965). 

[49] A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev, Integrals and Series (Gordon &; Breach, 
New York, 1989-1990), Vols. 1-3. 

15 



